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Abstract Particles with resonant short-range interactions have universal properties that do 
not depend on the details of their structure or their interactions at short distances. In the 
three-body system, these properties include the existence of a geometric spectrum of three-body 
Efimov states and a discrete scaling symmetry leading to log-periodic dependence of observables 
on the scattering length. Similar universal properties appear in the four-body system and 
possibly higher-body systems as well. For example, universal four-body states have recently 
been predicted and observed in experiment. These phenomena are often referred to as "Efimov 
Physics" . We review their theoretical description and discuss applications in different areas of 
physics with a special emphasis on nuclear and particle physics. 
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The scattering of particles with sufficiently low kinetic energy is determined by 
their S-wave scattering length a. This is the case if their de Broglie wavelengths 
are large compared to the range of the interaction. 

Generically, the scattering length a is comparable in magnitude to the range 
I of the interaction: \a\ ~ I. In exceptional cases, the scattering length can 
be much larger in magnitude than the range: \a\ 3> t. Such a large scattering 
length requires the fine-tuning of a parameter characterizing the interactions to 
the neighborhood of a critical value at which a diverges to ±00. If the scattering 
length is large, the particles exhibit properties that depend on a but are insensitive 
to the range and other details of the short-range interaction. These properties are 
universal in the sense that they apply equally well to any nonrelativistic particle 
with short range interactions that produce a large scattering length [TJ[2] . 

For example, in the case of equal-mass particles with mass m and a > 0, there 
is a two-body bound state near the scattering threshold with binding energy 
Bd = ft 2 1 (ma 2 ). The corrections to this formula are suppressed by powers of I /a. 
This bound state corresponds to a pole of the two-particle scattering amplitude 
at E = —Bd- If the scattering length is negative, there is a universal virtual 
state which corresponds to a pole on the unphysical second sheet in the complex 
energy plane. 

The key evidence for universal behavior in the three-body system was the 
discovery of the Efimov effect in 1970 13]. In the unitary limit l/a — I 0, the two- 
body bound state is exactly at the two-body scattering threshold E = 0. Efimov 
showed that in this limit there are infinitely many, arbitrarily-shallow three-body 

(n) 

bound states whose binding energies B\ have an accumulation point at E = 0. 
The Efimov effect is just one aspect of universal properties in the three-body 
system. It has universal properties not only in the unitary limit, but whenever 
the scattering length is large compared to the range I. In particular, the log- 
periodic scattering length dependence of observables is a unique consequence of 
Efimov physics. 

Although well established theoretically, the umambiguous identification of Efi- 
mov states in nature is difficult since typical systems are neither in the unitary 
limit nor can the scattering length be varied. The probably simplest example in 
nuclear physics is the triton. The triton can be interpreted as the ground state of 
an Efimov spectrum in the pnn-system with total spin J = 1/2. Since the ratio 
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£/a is only about one third, the whole spectrum contains only one state but the 
low-energy properties of the triton can be described in this scenario. A promising 
system for observing several Efimov states is 4 He atoms, which have a scattering 
length that is more than a factor of 10 larger than the range of the interaction. 
Calculations using accurate potential models indicate that the system of three 
4 He atoms has two three-body bound states or trimers. The ground-state trimer 
can be interpreted as an Efimov state, and it has been observed in experiments 
involving the scattering of cold jets of 4 He atoms from a diffraction grating |4j. 
The excited trimer is universally believed to be an Efimov state, but it has not 
yet been observed. 

The rapid development of the field of cold atom physics has opened up new 
opportunities for the experimental study of Efimov physics. This is made possible 
by two separate technological developments. One is the technology for cooling 
atoms to the extremely low temperatures where Efimov physics plays a crucial 
role. The other is the technology for controlling the interactions between atoms. 
By tuning the magnetic field to a Feshbach resonance, the scattering lengths of 
the atoms can be controlled experimentally and made arbitrarily large. Both 
developments were crucial in recent experiments that provided the first indirect 
evidence for the existence of Efimov states in ultracold atoms [5j. 

Overviews of Efimov physics in ultracold atomic gases can be found in Refs. [TJ 
[2J6] . In this review, we focus on universal aspects and Efimov states in nuclear and 
particle physics. Even though the scattering length can not be varied, there are 
many systems close to the unitary limit where Efimov physics is relevant. They 
include few-nucleon systems, halo nuclei, and weakly bound hadronic molecules. 
These systems can be described in a universal effective field theory (EFT) that 
implements an expansion around the unitary limit. Three-body bound states can 
be interpreted as Efimov trimers. 

In the next section we will review the physics of the Efimov effect starting 
with a brief account of the history. In the following sections, we will discuss 
applications in nuclear and particle physics. We will end with a summary and 
outlook. 
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2.1 History 

The first hints on universal behavior in the three-body system came from the 
discovery of the Thomas collapse in 1935 [7] which is closely related to the Efimov 
effect. Thomas studied the zero-range limit for potentials with a single two-body 
bound state with fixed energy. Using a variational argument, he showed that the 
binding energy of the deepest three-body bound state diverges to infinity 
in this limit. Thus the spectrum of three-body bound states is unbounded from 
below. 

Further progress occured by applying the zero-range limit to the three-nucleon 
system. An integral equation for S-wave neutron-deuteron scattering using zero- 
range interactions was derived by Skorniakov and Ter-Martirosian in 1957 (§]. 
For the spin-quartet channel, this integral equation has no bound state solutions 
and is well behaved. In the spin-doublet channel, however, it has solutions for 
arbitrary energy [§J, including bound state solutions. If the solution is fixed 
by requiring some three-body energy, the resulting equation still has a discrete 
spectrum that extends to minus infinity in agreement with the earlier result of 
Thomas 10, IT]. Although a prediction for the spin-doublet neutron-deuteron 
scattering length was obtained using the triton binding energy as input |12] , 
most work afterwards focused on finite range forces which avoid this pathology 
at high energies. 

In 1970, Efimov realized that one should focus on the physics at low energies, 
E <C ft 2 /(ml 2 ), and not on the deepest states. In this limit, where zero-range 
forces are adequate, he found some surprising results [3j. He pointed out that 
when \a\ is sufficiently large compared to the range I of the potential, there is a 
sequence of three-body bound states whose binding energies are spaced roughly 
geometrically in the interval between fi 2 /(ml 2 ) and Ti 2 /(ma 2 ). As \a\ is increased, 
new bound states appear in the spectrum at critical values of a that differ by 
multiplicative factors of e n ^ s °, where so depends on the statistics and the mass 
ratios of the particles. In the case of spin-doublet neutron-deuteron scattering 
and for three identical bosons, sq is the solution to the transcendental equation 

7TS 8 . 7TS - . 

sn cosh = — — smn . 1 

2 ^3 6 V ; 

Its numerical value is so ~ 1.00624, so e n ^ s ° ~ 22.7. As \a\/l —> oo, the asymp- 



6 Hammer, Platter 

totic number of three-body bound states is 

*T— ?b|H. (2) 

In the limit a — > ±00, there are infinitely many three-body bound states with 
an accumulation point at the three-body scattering threshold with a geometric 
spectrum: 

B (n) = ( e -^/^n-n, h 2 K 2j m ^ ^ 

where m is the mass of the particles and is the binding wavenumber of the 
branch of Efimov states labeled by n*. The geometric spectrum in ^ is the 
signature of a discrete scaling symmetry with scaling factor e n ^ s ° ~ 22.7. It is 
independent of the mass or structure of the identical particles and independent 
of the form of their short-range interactions. The Efimov effect can also occur in 
other three-body systems if at least two of the three pairs have a large S-wave 
scattering length but the numerical value of the asymptotic ratio may differ from 
the value 22.7. 

A formal proof of the Efimov effect was subsequently given by Amado and 
Noble 13 141. The Thomas and Efimov effects are closely related. The deepest 



three-body bound states found by Thomas's variational calculation can be iden- 
tified with the deepest Efimov states j To] . The mathematical connection of the 



Efimov effect to a limit cycle was discussed in 16 . 

The universal properties in the three-body system with large scattering length 
are not restricted to the Efimov effect. The dependence of three-body observ- 
ables on the scattering length or the energy is characterized by scaling behavior 



modulo coefficients that are log-periodic functions of a 17,18 . This behavior 
is characteristic of a system with a discrete scaling symmetry. We will refer to 
universal aspects associated with a discrete scaling symmetry as Efimov physics. 

In 1981, Efimov proposed a new approach to the low-energy few-nucleon prob- 
lem in nuclear physics that, in modern language, was based on perturbation 
theory around the unitary limit |19| . Remarkably, this program works reason- 
ably well in the three- nucleon system at momenta small compared to M n . The 
Efimov effect makes it necessary to impose a boundary condition on the wave 
function at short distances. The boundary condition can be fixed by using either 
the spin-doublet neutron-deuteron scattering length or the triton binding energy 
as input. If the deuteron binding energy and the spin-singlet nucleon- nucleon 
scattering length are used as the two-body input and if the boundary condition 
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is fixed by using the spin-doublet neutron-deuteron scattering length as input, 
the triton binding energy is predicted with an accuracy of 6%. The accuracy of 
the predictions can be further improved by taking into account the effective range 



as a first-order perturbation 20 . Thus the triton can be identified as an Efimov 
state associated with the deuteron and the spin-singlet virtual state being a pn 



state with large scattering length 19 



In the three-nucleon system, this program was implemented within an effec- 



tive field theory framework by Bedaque, Hammer, and van Kolck 21-23 . In 
Ref. [23], they found that the renormalization of the effective field theory re- 
quires a S'C/(4)-symmetric three-body interaction with an ultraviolet limit cycle. 
The three-body force depends on a parameter A* that is determined through a 
renormalization condition that plays the same role as Efimov's boundary condi- 
tion. 5{7(4)-symmetry was introduced by Wigner in 1937 as generalization of the 
SU(2) x SU(2) spin-isospin symmetry, allowing for a mixing of spin and isopin 
degrees of freedom in symmetry transformations [24] . It is satisfied to a high 
degree in the energy spectra of atomic nuclei. Exact Wigner symmetry requires 
the S-wave scattering lengths in the spin-triplet and spin-singlet channels to be 
equal. However, if the two-body scattering lengths are large, it is a very good 
approximation even if they are different since the symmetry-breaking terms are 



proportional to the inverse scattering lengths 25 



This effective field theory (EFT) is ideally suited to calculating corrections to 
the universal results in the scaling limit. Its application to few-nucleon physics 
will be discussed in Section 03 



2.2 Hyperspherical Methods 

Coordinate space methods have proven to be a valuable tool for the analysis of 



the three-body problem with short-range interactions 17 26 27J . In this subsec- 
tion we will introduce the hyperspherical approach that has been used to obtain 
important results about the Efimov spectrum. The material of this subsection 
is based on the discussion of the hyperradial formalism in the review of Ref. [l]. 
For three particles of equal mass the Jacobi coordinates are defined as: 



T k - -(Ti+Tj) , (4) 
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where the triple (ijk) is a cyclic permutation of the particle indices (123). The 
hyperradius R and hyperangle a^ are then defined by 

1 



R 



1 2 i 2 2 



afc = arctan 



V3\rij\ 



2r 



Mil 



(5) 



In the center-of-mass system the Schrodinger equation in hyperspherical coordi- 
nates is given by 

/ A? , 



with 



Tr + T a . + 



Tr 
A 2 

Mi 



Mi 
2mi? 2 



+ V(12, «) a, 0) = a, fi), 



iT 5/2 , 
2m V <9i? 2 4i? 2 



9 15 

+ 



i? 5 / 2 , 



h 2 



l 



(9 2 



2mR 2 sin 2a \ <9a 2 

T 2 T 2 

, ^Mi 



4 sin 2a, 



sin 2 afc cos 2 a^ 



(6) 

(7) 
(8) 
(9) 



where ft = (Oij , $ij , Ok,ij i <t>k,ij) an d the Ls that appear in Eq. ([9j) are the usual 
angular-momentum operators with respect to these angles. 

We will assume that the potential V depends only on the magnitude of the 
inter-particle separation and write 



F(n, r 2 , r 3 ) = V(r 12 ) + V(r 23 ) + V(r 31 ). 



(10) 



We now employ the usual Faddeev decomposition of ip for three identical bosons 
and neglect subsystem angular momentum 



V(R, a, O) = tp(R, ai) + ip(R, a 2 ) + ip(R, a 3 ). 



(11) 



The solution of the corresponding Faddeev equation can then be expanded in 
a set of eigenfunctions of the hyperangular operator, i.e. 

1 



ip(R, a) 



R 5 / 2 sin(2a) 



a 



(12) 



This leads to separate differential equations for the hyperangular functions <j) n 
and the hyperradial functions f n . We obtain in particular for the hyperradial 
functions 



EUR) 



It 
2m 



+ 



15 



dR 2 AR 2 



V n (R) 




UR) 



2Pnm{R)-Q^, + Qnm{R) 



fm(R), (13) 
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with the hyperradial potential V n {R) defined by 

V n (R) = (X n {R)-A) 



ti 2 



(14) 



2mR 2 ' 

and P nm (R) and Q n m(R) potentials that induce coupling between different hy- 
perradial channels [I]. 

For hyperradii R which are much larger than the range £ over which V is non- 
zero, the solution of the equation for the hyperangular function (j) n for large a 
is 



^ high) («) 



sin 



A„ ( | - a 



(15) 



On the other hand, for R 3> £ and small a, it can be shown that the solution 
for the hyperangular part can be written as 



8a 



n low) (a) = A^ Q (V2Ra) - T sm 

V3 



*-§) 



(16) 



where A is a constant and ■i/'o is the zero-energy solution to a two-body Schroedinger 
equation with the two-body potential V 

sin(kr + 5(k)) sm5(k) 



ipk(r) 

As k — > this yields ipo(r) 



A: 



[cos(fcr) + cot 5 sin(fer)] . 



(17) 



a, and we can use this asymptotic two-body 



wave function in Eq. (16). This gives 

4 low) (a) = A(V2Ra-a) 



8a 

71 



sin 



An 



(18) 



But, since V = in this region, this result must be consistent with Eq. (15). 
This is achieved by the choice 



A 



1 



sin 



7T 



(19) 



which ensures that (j) n {a) is continuous across the boundary between "low" and 
"high" solutions at a ~ £/R, and the condition 



cos 



An" 



-7T 



: sin 



7T 



sin 



7T\ R 



a 



(20) 



y/3K V v "6, 

on A n , which ensures that cp n (a) has a continuous first derivative as a — > l/R. 
We note that if these equations are satisfied A n , and hence (f> n , is independent of 



R for R <^ \a\. Indeed, as long as Eqs. (20) and (19) are satisfied the form (15) is 



the result for <fi for all a such that a > l/R. Solving Eq. (20) in the limit R <C |o| 
we find the lowest eigenvalue 



An 



-s(iW) 



(21) 
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with so = 1.00624.... This is the only negative eigenvalue, and hence only this 
channel potential is attractive. So, if we now focus on the unitary limit, where 
\a\ — > oo, we have Aq = — Sq. Since it can also be shown that the coupling 



potentials P nm and Q n m vanish in this regime, the hyperradial equation ( 13 ) in 
the lowest channel becomes 



h 2 ( d 2 jj+j 

2m \ dR 2 R 2 



fo(R) = Efo(R). (22) 



This equation will hold for R S> i- If we desire a solution for negative E the 
requirement of normalizability for /o mandates that 

f(°\R) = VRK iS0 (V2KR), (23) 

where the superscript (0) indicates that we are working in the unitary limit, while 
the subscript refers to the solution for the hyperchannel corresponding to Ao, 
which is the only one that supports bound states. The binding energy of these 



bound states is related to the k of Eq. ( 23 ) by 



B t = . 24 

m 

Since the attractive 1/R 2 potential produces a spectrum that is unbounded from 
below some other short-distance physics is needed in order to stabilize the system. 
If the two-body potential is known this short-distance physics is provided by the 
two-body potential V, that becomes operative for R ~ I. But an alternative 



approach is to add an additional term to Eq. (22) that summarizes the impact of 
the two-body V. Here we take this potential to be a surface delta function at a 
radius 1/A [28] 

V SR (R) = H (A)A 2 5 \R- , (25) 

with Hq adjusted as a function of A such that the binding energy of a particular 
state, say B^ 1 *^ (with a corresponding k*, given by (24)), is reproduced. Note 



that since Vsr is operative only at small hyperradii R ~ 1/A it corresponds to 



a three-body force. (See Ref. 29 for a realization of this in a momentum-space 
formalism.) In physical terms we anticipate t ~ 1/A, since we know that once we 
consider hyperradii of order 1 /A the potential V starts to affect the solutions. 

Given that our focus is on predictions of the theory that are independent of 
details of V we can consider the extreme case and take the limit £ — > 0. In this 
limit the form of Ki SQ as R — > guarantees that once Hq is fixed to give a bound 
state at B\ , the other binding energies in this hyperradial eigenchannel form 
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a geometric spectrum. Namely, B^ = h 2 K^/m with 



't 

n—n* 



(26) 



with ra* the index of the bound state corresponding to «;*. Eq. (26) will hold for 
all K n such that n n -C A. (Note that now the subscript on k denotes the index 
of the bound state in adiabatic channel zero.) The continuous scale invariance 
of the 1/R 2 potential has been broken down to a discrete scale invariance by 
the imposition of particular short-distance physics on the problem through the 



short-distance potential (25) |28|. 

It was subsequently shown that the discrete scale invariance of the three-body 
wavefunction which is exact in the limit of infinite scattering length and zero range 
also has implications for finite range. A perturbative calculation of the effect of 
a finite effective range on the bound state spectrum showed that the spectrum 



remains unchanged 30 and corrections to binding energies are of order (ro / a) 



2 



where ro is the effective range of the interaction. 



2.3 Efimov Spectrum 

The hyperspherical methods discussed above can also be used to obtain the bind- 
ing energy spectrum at finite scattering length. The short-distance boundary con- 
dition that was used to fix the binding energy in the unitary limit also determines 
the bound state spectrum at finite scattering length. The binding momentum 
introduced above can therefore be considered as a convenient parameter that 
determines the value of all universal few-body observable of the corresponding 
universality class. 

The exact discrete scaling symmetry observed in the limit of infinite scattering 
length also exists if is kept fixed and a and other variables such as the energy 
are rescaled 

— y , a^SZ 1 a , S^ 2m E . (27) 

Observables such as binding energy and cross sections scale with integer powers 
of Sq = exp(7r/so) under this symmetry. For example, the binding energy of an 
Efimov trimer which is a function of a and k* scales 

B[ n) (5 m a, «») = S^ 2m B[ n - m) (a, «,) . (28) 

This implies for positive scattering length 

B[ n \a,^) = F n (2s ln(arc*))^ . (29) 
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Figure 1: The Efimov plot for the three-body problem. We show K = 
sgn(E)(m\E\) 1 ^ 2 /h versus the inverse scattering length. The allowed regions for 
three-atom scattering states and atom-dimer scattering states are labeled AAA 
and AD, respectively. The heavy lines labeled T are two of the infinitely many 
branches of Efimov states. The cross-hatching indicates the threshold for scat- 
tering states. States along the vertical dashed line have a fixed scattering length. 



The function F n parametrizes the scattering length dependence of all Efimov 
trimers exactly in the limit of vanishing range. The function F n satisfies 

F n {x + 2mvr) = {e~^ ' sa ) m F n . m (x) . (30) 

The scattering length dependence of the bound state spectrum is shown in Fig. 
[lj We plot the quantity K = sgri{E){m\E\) 1 / 2 /h against the inverse scattering 
length. For bound states K corresponds to the binding momentum. The lines 
denote Efimov trimers below the threshold. The threshold for scattering states is 
denoted by the hatched area. Only a few of the infinitely many Efimov branches 
are shown. A given physical system has a fixed scattering length value and is 
denoted by the vertical dashed line. Changing the parameter by a factor 
Sq corresponds to multiplying each branch of trimers with this factor without 
changing their shapes. One important result is that three-bound states exist for 
positive and negative scattering length. This is remarkable for the latter case 
since the two-body subsystem is unbound for a < 0. At a negative scattering 
length that we denote with a^, a bound state with a give has zero binding 
energy. As the scattering length is increased the binding energy gets larger until 
it crosses the atom-dimer threshold at the positive scattering length a*. The 



Efimov States in Nuclear and Particle Physics 13 

quantities a* and a* can also be used to quantify a universality class of Efimov 
states. 



2.4 Universal Properties 

Other calculable observables will also display the discrete scaling symmetry that 
we just discussed for the bound state spectrum. The atom-dimer cross section 
fulfils for example the constraint 

a AD (S^ 2m E; 5 m a, k«) = S 2m a AD (E; a, fc*) (31) 

under rescaling. At E = the cross section is related to the atom-dimer scattering 
length (Tad = 47r|aAD| 2 - This implies that the atom-dimer scattering length can 
be written as 

oad = f(2s ln(an*))a , (32) 

where f(x) is periodic function with period 2tt. 

An observable that has been crucial for the experimental detection of Efimov 
physics in ultracold atoms is the three-body recombination rate. In ultracold 
gases, atoms can undergo inelastic three-body collisions in which a two-body 
bound state is formed. The dimer and remaining atom gain kinetic energy in this 
process and can leave the atomic trap. Such processes lead to a measurable loss 
of particles in the atomic trap. For negative scattering length, atoms can only 
recombine into deep dimers that have binding energy of order h 2 /(mR 2 ). For 
positive scattering length, atoms can recombine into shallow dimers with binding 
energy Ti 2 /{ma 2 ) and deep dimers. The recombination rate that is a measure for 
the loss rate of atoms scales as Ha 4 /m times a log-periodic coefficient that is a 
function of a and . The recombination rate constant a can therefore be written 
as 

a (a) = <?(2so ln(aK*)) . (33) 

m 

The analytic form of the function g(x) is known and can be found in [TJ[6]. 
Here we want to focus on the qualitative features as shown in Fig. [2] for positive 
scattering length. At positive scattering length interference effects lead to log- 
periodically spaced minima in the recombination rate. At negative scattering 
length free atoms can only recombine into deep dimers. This process will be 
enhanced dramatically whenever an Efimov trimer is at threshold. 

A different perspective on Efimov physics can be gained by keeping the two- 
body scattering length fixed and varying the three-body parameter. As a result 
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all three-body observables are correlated and will lead to correlation lines when 
plotted against each other. One of them (well-known from nuclear physics) is 
the Phillips line [3l], the correlation that results when the trimer binding en- 
ergy is plotted against the atom-dimer scattering length. Such correlations had 
been observed frequently in nuclear physics, where different phase shift equivalent 
potentials were employed in few-body calculations. 

Since one three-body parameter is required for a description of the three-body 
system with zero-range interactions it is natural to ask how many parameters are 
needed for calculations in the n-body system. A first step towards answering this 
question was performed in Ref. [32]. The authors of this work showed that the 
two-body scattering length and one three-body parameter are sufficient to make 
predictions for four-body observables. Results in a more detailed analysis [33] also 
lead to the conclusion that every trimer state is tied to two universal tetramer 
states with binding energies related to the binding energy of the next shallower 
trimer. In the unitary limit 1/a = 0, the relation between the binding energies 
was found as: 

Bf ] ^5B t and rj 1.01 J3 t . (34) 

where B^ denotes the binding energy of the deeper of the two tetramer states 
and B^ the shallower of the two. 

A recent calculation by von Stecher, d'Incao and Greene [34] supports these 
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findings and extends them to higher numerical accuracy. For the relation between 
universal three- and four-body bound states in the unitary limit, they found 

Bf ] ^4.57B t and « 1.01 B t , (35) 



which is consistent with the results given in Eq. ( 34 ) within the numerical accu- 
racy. 

The results obtained by the Hammer and Platter in Ref. [33] were furthermore 
presented in the form of an extended Efimov plot, shown in Fig. [3j Four-body 
states must have a binding energy larger than the one of the deepest trimer 
state. The corresponding threshold is denoted by lower solid line in Fig. [3j The 
threshold for decay into the shallowest trimer state and an atom is indicated 
by the upper solid line. At positive scattering length, there are also scattering 
thresholds for scattering of two dimers and scattering of a dimer and two particles 
indicated by the dash-dotted and dashed lines, respectively. The vertical dotted 
line denotes infinite scattering length. 

An extended version of this four-body Efimov plot was also presented by von 
Steelier, d'Incao and Greene in Ref. |34| . They calculated more states with higher 
numerical accuracy and extended the calculation of the four-body states to the 
thresholds where they become unstable. From these results they extracted the 
negative values of the scattering lengths at which the binding energies of the 
tetramer states become zero and found 

a\ Q « 0.43< and a\ x « 0.92< . (36) 

These numbers uniquely specify the relative position of three- and four-body 
recombination resonances. This was the key information for the subsequent ob- 
servation of these states in ultracold atoms by Ferlaino et al. |35|. 

Calculations for larger number of particles using a model that incorporates the 
universal behavior of the three-body system were carried our by von Stecher [36] . 
His findings indicate that there is at least one A-body state tied to each Efimov 
trimer and numerical evidence was also found for a second excited 5-body state. 



2.5 Observation in Ultracold Atoms 

The first experimental evidence for Efimov physics in ultracold atoms was pre- 
sented by Kraemer et al. (5 in 2006. This group used 133 Cs atoms in the lowest 
hyperfine spin state. They observed resonant enhancement of the loss of atoms 
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Figure 3: The extended Efimov plot for the four-body problem. We show K = 
sgn(E)(m\E\) 1 / 2 /fi versus the inverse scattering length. Both quantities are given 
in arbitrary units. The circles and triangles indicate the four-body ground and 
excited states, respectively, while the lower (upper) solid lines give the thresholds 
for decay into a ground state (excited state) trimer and a particle. The dash- 
dotted (dashed) lines give the thresholds for decay into two dimers (a dimer and 
two particles). The vertical dotted line indicates infinite scattering length. 



from three-body recombination that can be attributed to an Efimov trimer cross- 
ing the three-atom threshold. Kraemer et al. also observed a minimum in the 
three-body recombination rate that can be interpreted as an interference effect 
associated with Efimov physics. In a subsequent experiment with a mixture of 
133 Cs atoms and dimers, Knoop et al. observed a resonant enhancement in the 
loss of atoms and dimers [37]. This loss feature can be explained by an Efimov 



trimer crossing the atom-dimer threshold 38 . The most exciting recent devel- 
opments in the field of Efimov physics involve universal tetramer states. Ferlaino 
et al. observed two tetramers in an ultracold gas of 133 Cs atoms [35] that con- 



firm the results by Hammer and Platter [32,33 and von Stecher, D'Incao and 
Greene [34] . 

Recent experiments with other bosonic atoms have provided even stronger 
evidence of Efimov physics in the three- and four-body sectors. Zaccanti et al. 
measured the three-body recombination rate and the atom-dimer loss rate in a 
ultracold gas of 39 K atoms 39 . They observed two atom-dimer loss resonances 



and two minima in the three-body recombination rate at large positive values 
of the scattering length. The positions of the loss features are consistent with 
the universal predictions with discrete scaling factor 22.7. They also observed 
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loss features at large negative scattering lengths. Barontini et al. obtained the 
first evidence of the Efimov effect in a heteronuclear mixture of 41 K and Rb 
atoms 



40 . They observed 3-atom loss resonances at large negative scattering 



lengths in both the K-Rb-Rb and K-K-Rb channels, for which the discrete scaling 
factors are 131 and 3.51 x 10 5 , respectively. Gross et al. measured the three- 



body recombination rate in an ultracold system of 7 Li atoms 41 . They observed 
a 3-atom loss resonance at a large negative scattering length and a three-body 
recombination minimum at a large positive scattering length. The positions of 
the loss features, which are in the same universal region on different sides of a 
Feshbach resonance, are consistent with the universal predictions with discrete 
scaling factor 22.7. Pollack et al. measured the three-body recombination in a 
system of 7 Li atoms in a hyperfine state different from the system considered by 



Gross et al. 42 . They observed a total of 11 three- and four-body loss features. 
The features obey the universal relations on each side of the Feshbach resonance 
separately, however, a systematic error of 50 % is found when features on different 
sides of the Feshbach resonance are compared. 

Efimov physics has also been observed in three-component systems of 6 Li atom. 
For the three lowest hyperfine states of 6 Li atoms, the three pair scattering lengths 
approach a common large negative value at large magnetic fields and all three 
have nearby Feshbach resonances at lower fields that can be used to vary the 
scattering lengths [43| . The first experimental studies of many-body systems of 
6 Li atoms in the three lowest hyperfine states have recently been carried out 
by Ottenstein et al. [44] and by Huckans et al. [45]. Their measurements of 
the three-body recombination rate revealed a narrow loss feature and a broad 
loss feature in a region of low magnetic field. Theoretical calculations of the 
three-body recombination rate supported the interpretation that the narrow loss 



feature arises from an Efimov trimer crossing the 3-atom threshold [46 - 48 . Very 
recently, another narrow loss feature was discovered in a much higher region of 
the magnetic field by Williams et al. [49| and by Jochim and coworkers. Williams 
et al. used measurements of the three-body recombination rate in this region to 
determine the complex three-body parameter that governs Efimov physics in this 
system. This parameter, together with the three scattering lengths as functions 
of the magnetic field, determine the universal predictions for 6 Li atoms in this 
region of the magnetic field. 
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The properties of hadrons and nuclei are determined by quantum chromodynam- 
ics (QCD), a non-abelian gauge theory formulated in terms of quark and gluon 
degrees of freedom. At low energies, however, the appropriate degrees of free- 
dom are the hadrons. Efimov physics and the unitary limit can serve as a useful 
starting point for effective field theories (EFTs) describing hadrons and nuclei 
at very low energies. For convenience, we will now work in natural units where 

h = c = i. 

In nuclear physics, there are a number of EFTs which are all useful for a cer- 



tain range of systems 50-52 . At very low energies, where Efimov physics plays 
a role, all interactions can be considered short-range and even the pions can be 
integrated out. This so-called "pionless EFT" is formulated in an expansion of 
the low- momentum scale M\ ow over the high- momentum scale M^igh- It can be 
understood as an expansion around the limit of infinite scattering length or equiv- 
alently around threshold bound states. Its breakdown scale is set by one-pion 
exchange, M^g^ ~ M w , while Mi ow ~ 1/a ~ k. For momenta k of the order of 
the pion mass M n , pion exchange becomes a long-range interaction and has to be 
treated explicitly. This leads to the chiral EFT whose breakdown scale M^gh is 
set by the chiral symmetry breaking scale A x . The pionless theory relies only on 
the large scattering length and is independent of the short-distance mechanism 
generating it. This theory is therefore ideally suited to unravel universal phe- 



nomena driven by the large scattering length such as limit cycle physics 153 54 



and the Efimov effect [3j. In this section, we will focus on the aspects of nuclear 
effective field theories related to Efimov physics. For more complete overviews 
of the application of effective field theories to nuclear phenomena in general, a 



number of excellent reviews are available |50f 52 , 55 



3.1 Few-Nucleon System 



In the two-nucleon system, the pionless theory reproduces the well known effective 
range expansion in the large scattering length limit. The renormalized S-wave 
scattering amplitude to next-to-leading order in a given channel takes the form 



T 2 (k) 



4tt 



1 



1 



r k 2 /2 
-1/a — ik 



+ 



(37) 



m —1/a — ik 

where k is the relative momentum of the nucleons and the dots indicate correc- 
tions of order {Mi ow j 'Mhigh) 2 for typical momenta k ~ M\ ow . In the language 
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of the renormalization group, this corresponds to an expansion around the non- 
trivial fixed point for 1/a = [56 , 57 . The pionless EFT becomes very useful 



in the two-nucleon sector when external currents are considered and has been 
applied to a variety of electroweak processes. These calculations are reviewed in 
detail in Refs. |50[[5l). 

Here we focus on the three-nucleon system. It is convenient (but not manda- 
tory) to write the theory using so-called "dimeron" auxiliary fields 1 58 1 . We need 
two dimeron fields, one for each S-wave channel: (i) a field U with spin (isospin) 
1 (0) representing two nucleons interacting in the 3 Si channel (the deuteron) and 
(ii) a field s a with spin (isospin) (1) representing two nucleons interacting in 


4m 
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(38) 



where i,j are spin and a, b are isospin indices while gt, g s , At, A s and G3 are the 
bare coupling constants. The Pauli matrices Oi (r a ) operate in spin (isospin) 
space, respectively. This Lagrangian goes beyond leading order and already 
includes the effective range terms. The coupling constants gt, At, g s , A s are 
matched to the scattering lengths a a and effective ranges r^ a in the two channels 
(a = s,t). Alternatively, one can match to the position of the bound state/virtual 
state pole 7^ in the T-matrix instead of the scattering length which often im- 
proves convergence 



59 



The term proportional to G3 constitutes a Wigner-5'C/(4) symmetric three- 
body interaction. It only contributes in the spin-doublet S-wave channel. When 
the auxiliary dimeron fields t% and s a are integrated out, an equivalent form 
containing only nucleon fields is obtained. At leading order when the effective 
range corrections are neglected, the spatial and time derivatives acting on the 
dimeron fields are omitted and the field is static. The coupling constants g a and 
A a , a = s,t are then not independent and only the combination g\j A a enters 
in observables. This combination can then be matched to the scattering length 
or pole position. 

The simplest three-body process to consider is neutron-deuteron scattering 
below the breakup threshold. In order to focus on the main aspects of renor- 
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Figure 4: The integral equation for the particle-dimeron scattering amplitude. 
The single (double) line indicates the particle (dimeron) propagator. 

malization, we suppress all spin-isospin indices and complications from coupled 
channels in the three-nucleon problem. This is equivalent to a system of three 
spinless bosons with large scattering length. If the scattering length is posi- 
tive, the particles form a two-body bound state analog to the deuteron which we 
generically call dimeron. The leading order integral equation for particle-dimeron 
scattering is shown schematically in Fig. [4j Projected on total orbital angular 
momentum L = 0, it takes the following form: 

16 4 f A 

T 3 (k,p;E) = M(k,p;E) + - dqq 2 T 3 (k,q; E) 
3a 7r J 

M(q,p; E) 



= , (39) 

-1/a + y / 3q 2 /A-mE -ie 

where the inhomogeneous term reads 

1 , fk 2 + kp + p 2 -mE\ H(A) lA , 

Here, if (A) is a running coupling constant that determines the strength of the 
three-body force G 3 (A) = 2mH(A)/A 2 and A is a UV cutoff introduced to reg- 
ularize the integral equation. Note that the three-body force is enhanced and 
enters already at leading order in this theory. The magnitude of the incom- 
ing (outgoing) relative momenta is k (p) and E = 3k 2 /(4m) — l/(ma 2 ). The 
on-shell point corresponds to k = p and the phase shift can be obtained via 



kcotd = 1/T 3 (k,k; E) + ik. For if = and A ->• oo, Eq. (|39j) reduces to the STM 
equation of Skorniakov and Ter-Martirosian [8J. It is well known that the STM 
equation has no unique solution [9] . The regularized STM equation has a unique 
solution for any given (finite) value of the ultraviolet cutoff A but the solution 
strongly depends on the value of A. In the EFT framework, cutoff independence 



of the amplitude is achieved by an appropriate "running" of if (A) 29 , 60 



= cos [s ln(A/A^ + arctan s ] 
cos[so h^A/A*) — arctanso] 
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where A* is a dimensionful three-body parameter generated by dimensional trans- 
mutation. The dependence of the three-body coupling H on the cutoff A is shown 



< 
x 





Figure 5: The three-body coupling H as a function of the cutoff A for a fixed value 
of the three-body parameter A*. The solid line shows the analytical expression 
(41), while the dots show results from the numerical solution of Eq. (39). 



in Fig. [5] for a fixed value of the three-body parameter A*. The solid line shows 



the analytical expression (41), while the dots show results from the numerical 



solution of Eq. (39). A good agreement between both methods is observed, in- 



dicating that the renormalization is well under control. Adjusting A* to a single 
three-body observable allows to determine all other low-energy properties of the 
three-body system* 1 Because H(A) in Eq. (41) vanishes for certain values of 



the cutoff A it is possible to eliminate the explicit three-body force from the 
equations by working with a fixed cutoff that encodes the dependence on A*. 
This justifies tuning the cutoff A in the STM equation to reproduce a three-body 
datum and using the same cutoff to calculate other observables as suggested by 
Kharchenko [61] . Equivalently, a subtraction can be performed in the integral 



equation 62 , 63 . In all cases one three-body input parameter is needed for the 



calculation of observables. 

The discrete scaling symmetry of the Efimov spectrum is manifest in the run- 
ning of the coupling H(A). The spectrum of three-body bound states of this EFT 
is exactly the Efimov spectrum. The integral equations for the three-nucleon 



problem derived from the Lagrangian ( 38 ) are a generalization of Eq. ( 39 ) . (For 



#1 Note that the choice of the three-body parameter A* is not unique, for alternative definitions 
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their explicit form and derivation, see e.g. Ref. |64|.) 

For S-wave nucleon-deuteron scattering in the spin-quartet channel only the 
spin-1 dimeron field contributes and the integral equation becomes [8,21, 22 



T^ /2 \p,k;E) = -^K(p,k)-- r dqq 2 D t ( q -E)K(p, q )T^ 2 \ q ,k;E) , 

11 Jo 

(42) 

where 

K( p k) = — In [ P 2 +P k + k2 - mE \ ( 43 N 
' pk \p 2 — pk + k 2 — mE J ' 

Dt(q;E) is the full spin-1 dimeron propagator and 74 ~ 45 MeV the deuteron 

pole momentum. This integral equation has a unique solution for A — > 00 and 

there is no three-body force in the first few orders. An S-wave three-body force 

is forbidden by the Pauli principle in this channel since all nucleon spins have to 

be aligned to obtain J = 3/2. The spin-quartet scattering phases /scotch 3 / 2 ) = 

l/Tf^'{k,k; E) + ik can therefore be predicted to high precision from two-body 

data alone. 

In the spin-doublet channel both dimeron fields as well as the three-body force 



in the Lagrangian ( 38 ) contribute 23 . This leads to a pair of coupled integral 
equations for the T-matrix. The renormalization of this equation is easily un- 
derstood in the unitary limit which corresponds to a Wigner SU(4) symmetry 
of the theory [24]. In the unitary limit these two integral equations decouple. 
One of the two equations has the same structure as the equation for the bosonic 



problem ( 39 ) , while the other one is similar to the equation in the quartet chan- 



nel (42). Thus, one needs one new parameter which is not determined in the 2N 
system in order to fix the (leading) low-energy behavior of the 3N system in this 
channel. This parameter corresponds to the SU (4) symmetric three-body force 



proportional to G3 in the Lagrangian (38) |23 . The three-body parameter gives 
a natural explanation of universal correlations between different three-body ob- 
servables such as the Phillips line, a correlation between the triton binding energy 
and the spin-doublet neutron-deuteron scattering length [3l] . These correlations 
are purely driven by the large scattering length independent of the mechanism 
responsible for it. If the spin-doublet neutron-deuteron scattering length is given, 
the triton binding energy is predicted. In this scenario the triton emerges as an 
Efimov state. The scenario can be tested by using the effective theory to predict 
other three-body observables. 

Higher-order corrections to the amplitude including the ones due to 2N effec- 
tive range terms can be included perturbatively. This was first done at NLO 
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Figure 6: Phase shifts for neutron-deuteron scattering below the deuteron 
breakup at LO (dash-dotted line), NLO (dashed line), and N 2 LO (solid line). 
The filled squares and circles give the results of a phase shift analysis and a 
calculation using AV18 and the Urbana IX three-body force, respectively. 



for the scattering length and triton binding energy in [20] and for the energy 



dependence of the phase shifts in [621. In Refs. 64,65 , it was demonstrated that 



it is convenient to iterate certain higher order range terms in order to extend the 
calculation to N 2 LO. Here, also a subleading three-body force was included as re- 
quired by dimensional analysis. More recently, Platter and Phillips showed using 
the subtractive renormalization that the leading three-body force is sufficient to 
achieve cutoff independence up to N 2 LO in the expansion in Mi^/M^igh 



66 



The 



results for the spin-doublet neutron-deuteron scattering phase shift at LO [23] , 
NLO 1 62 1, and N 2 LO |67| are shown in Fig. |6j There is excellent agreement with 
the available phase shift analysis and a calculation using a phenomenological NN 
interaction. From dimensional analysis, one would expect the subleading three- 
body force at N 2 LO. Whether there is a suppression of the subleading three-body 
force or simply a correlation between the leading and subleading contributions is 
not understood. 

Three-nucleon channels with higher orbital angular momentum are similar to 
the spin-quartet for S- waves and three-body forces do not appear until very high 



orders 68 . A general counting scheme for three-body forces based on the asymp- 



totic behavior of the solutions of the leading order STM equation was proposed 
by Griefihammer [69] . A complementary approach to the few-nucleon problem is 
given by the renormalization group where the power counting is determined from 
the scaling of operators under the renormalization group transformation |70) . 
This method leads to consistent results for the power counting 
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Figure 7: The Tjon line correlation as predicted by the pionless theory. The grey 
circles and triangles show various calculations using phenomenological potentials 
|81| . The squares show the results of chiral EFT at NLO for different cutoffs while 
the diamond gives the N 2 LO result (82}[83j. The cross shows the experimental 
point. 



Three-body calculations with external currents are still in their infancy. How- 
ever, a few exploratory calculations have been carried out. Universal properties of 



the triton charge form factor were investigated in Ref. 74 and neutron-deuteron 



radiative capture was calculated in Refs. 75-77 . Electromagnetic properties of 



the triton were recently investigated in Refs. 78 , 79 . This work opens the pos 



sibility to carry out accurate calculations of electroweak reactions at very low 
energies for astrophysical processes. 



The pionless approach has also been extended to the four-body sector 32 , 80 



In order to be able to apply the Yakubovsky equations, an equivalent effective 
quantum mechanics formulation was used. The study of the cutoff dependence of 
the four-body binding energies revealed that no four-body parameter is required 
for renormalization at leading order. As a consequence, there are universal corre- 
lations in the four-body sector which are also driven by the large scattering length. 
The best known example is the Tjon line: a correlation between the triton and 
alpha particle binding energies, Bt and B a , respectively. Of course, higher order 
corrections break the exact correlation and generate a band. In Fig. [7j we show 



this band together with some calculations using phenomenological potentials 81 



and a chiral EFT potential with explicit pions 82,83 . All calculations with in 



teractions that give a large scattering length must lie within the band. Different 
short-distance physics and/or cutoff dependence should only move the results 
along the band. This can for example be observed in the NLO results with the 
chiral potential indicated by the squares in Fig. [7] or in the few-body calculations 
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with the low- momentum NN potential Vi ow k carried out in Ref. 84 . The Vjowk 
potential is obtained from phenomenological NN interactions by integrating out 
high-momentum modes above a cutoff A but leaving two-body observables (such 
as the large scattering lengths) unchanged. The results of few-body calculations 
with Vi ow k are not independent of A but lie all close to the Tjon line (cf. Fig. 2 
in Ref. 

Another interesting development is the application of the Resonating Group 
Model to solve the pionless EFT for three- and four-nucleon systems [85]. This 
method allows for a straightforward inclusion of Coulomb effects. Kirscher et al. 
extended previous calculations in the four-nucleon system to next-to-leading order 
and showed that the Tjon line correlation persists. Moreover, they calculated the 
correlation between the singlet S-wave 3 He-neutron scattering length and the 
triton binding energy. Preliminary results for the halo nucleus 6 He have been 
reported in [86] . 

The pionless theory has also been applied within the no-core shell model ap- 
proach. Here the expansion in a truncated harmonic oscillator basis is used as the 
ultraviolet regulator of the EFT. The effective interaction is determined directly 
in the model space, where an exact diagonalization in a complete many-body ba- 



sis is performed. In Ref. 87 , the + excited state of He and the Li ground state 
were calculated using the deuteron, triton, alpha particle ground states as input. 
The first + excited state in 4 He is calculated within 10% of the experimental 
value, while the 6 Li ground state comes out at about 70% of the experimental 
value in agreement with the 30 % error expected for the leading order approxi- 
mation. These results are promising and should be improved if range corrections 
are included. Finally, the spectrum of trapped three- and four-fermion systems 



was calculated using the same method 88 . In this case the harmonic potential 



is physical and not simply used as an ultraviolet regulator. For an update on this 



work, see 89 



3.2 Quark Mass Dependence and Infrared Limit Cycle 

In the following, we discuss the possibility of an exact infrared limit cycle and the 
Efimov effect in a deformed version of QCD with quark masses slightly larger than 
their physical values. The quark mass dependence of the chiral NN interaction 
was calculated to next-to-leading order (NLO) in the chiral counting in Refs. |90[ 



91 . At this order, the quark mass dependence is synonymous to the pion mass 
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dependence because of the Gell-Mann-Oakes-Renner relation: M% = —(m u + 
m^) (0\uu\0) / F 2 , where (0|uu|0) ~ (—290 MeV) 3 is the quark condensate. In 
the following, we will therefore refer to the pion mass dependence instead of 
the quark mass dependence which is more convenient for our purpose. The pion 
mass dependence of the nucleon-nucleon scattering lengths in the 3 Si- 3 Di and 1 Sq 



channels as well as the deuteron binding energy were calculated in Refs. 90-92 



In principle, the pion mass dependence of the chiral NN potential is deter- 
mined uniquely. However, the extrapolation away from the physical pion mass 
generates errors. The dominating source are the constants C$t an d Ds,T which 
give the corrections to the LO contact terms oc M% and cannot be determined 
independently from fits to data at the physical pion mass. A smaller effect is 
due to the error in the LEC d±Q, which governs the pion mass dependence of qa- 
Both effects generate increasing uncertainties as one extrapolates away from the 
physical point. 



In the calculation of Ref. 91 , the size of the two constants D$ and Dt was 
constrained from naturalness arguments, assuming —3 < F%A 2 x D s ,t < 3 , where 
A x ~ 1 GeV is the chiral symmetry breaking scale. These bounds are in agree- 
ment with resonance saturation estimates 93 . The constant diQ was varied in the 



range d±Q = —0.91... — 1.76 GeV 2 94 . These ranges were used to estimate the 
extrapolation errors of two-nucleon observables like the deuteron binding energy 



and the spin-singlet and spin-triplet scattering lengths 91 . In the chiral limit the 
deuteron binding energy was found to be of natural size, B d ~ F 2 /m ~ 10 MeV. 
Note, however, that if larger uncertainties in the LECs D$ and Dt are assumed 
one cannot make a definite statement about the binding of the deuteron in the 
chiral limit |90[|92| . For pion masses above the physical value, however, all cal- 
culations show similar behavior. 

In Fig. [8l we show the inverse scattering lengths in the spin-triplet and spin- 



singlet channels from Ref. 91 together with some recent lattice results [97] . 
Figure [8] also shows that a scenario where both inverse scattering lengths vanish 
simultaneously at a critical pion mass of about 200 MeV is possible. For pion 
masses below the critical value, the spin-triplet scattering length is positive and 
the deuteron is bound. As the inverse spin-triplet scattering length decreases, the 
deuteron becomes more and more shallow and finally becomes unbound at the 
critical mass. Above the critical pion mass the deuteron exists as a shallow virtual 
state. In the spin-singlet channel, the situation is reversed: the "spin-singlet 
deuteron" is a virtual state below the critical pion mass and becomes bound 
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Figure 8: Inverse of the S-wave scattering lengths in the spin-triplet and spin- 
singlet nucleoli- nucleon channels as a function of the pion mass M n . Filled tri- 
angles and rectangles show the lattice calculations from Refs. 95, 96] and [97] , 
respectively. 



above. Based on this behavior Braaten ans Hammer conjectured that one should 
be able to reach the critical point by varying the up- and down-quark masses 
m u and rrid independently because the spin-triplet and spin-singlet channels have 



different isospin 54 . In this case, the triton would display the Efimov effect which 
corresponds to the occurrence of an infrared limit cycle in QCD. It is evident that 
a complete investigation of this issue requires the inclusion of isospin breaking 
corrections and therefore higher orders in the chiral EFT. However, the universal 
properties of the limit cycle have been investigated by considering specific values 
of Ds and Dt that lie within the naturalness bound and cause the spin-singlet 
and spin-triplet scattering lengths to become infinite at the same value of the 
pion mass. 

In Ref. |98| , the properties of the triton around the critical pion mass were 
studied for one particular solution with a critical pion mass M^ lt = 197.8577 
MeV. From the solution of the Faddeev equations, the binding energies of the 
triton and the first two excited states in the vicinity of the limit cycle were 
calculated for this scenario in chiral EFT. The binding energies are given in Fig. [9] 
by the circles (ground state), squares (first excited state), and diamonds (second 
excited state). The dashed lines indicate the neutron-deuteron (M n < M£ nt ) 
and neutron-spin-singlet-deuteron (M n > M£" ) thresholds where the three-body 
states become unstable. Directly at the critical mass, these thresholds coincide 
with the three-body threshold and the triton has infinitely many excited states. 
The solid lines are leading order calculations in the pionless theory using the pion 
mass dependence of the nucleon-nucleon scattering lengths and one triton state 
from chiral EFT as input. The chiral EFT results for the other triton states in 
the critical region are reproduced very well. The binding energy of the triton 
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Figure 9: Binding energies Bt of the triton ground and first two excited states as 
a function of M n . The circles, squares, and diamonds give the chiral EFT result, 
while the solid lines are calculations in the pionless theory. The vertical dotted 
line indicates the critical pion mass M£ nt and the dashed lines are the bound 
state thresholds. 
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Figure 10: Doublet neutron-deuteron scattering length a ' . in the critical region 
computed in the pionless EFT. The solid line gives the LO result, while the 
crosses and circles show the NLO and N 2 LO results. The dotted lines indicate 

1/2 

the pion masses at which a , diverges. 



ground state varies only weakly over the whole range of pion masses and is about 
one half of the physical value at the critical point. The excited states are strongly 
influenced by the thresholds and vary much more strongly. 

These studies were extended to N 2 LO in the pionless EFT and neutron-deuteron 
scattering observables in Ref. [99] . It was demonstrated that the higher order cor- 
rections in the vicinity of the critical pion mass are small. This is illustrated in 
Fig. 



10 



1/2 

where we show the doublet scattering length a_ . in the critical region. 
The solid line gives the LO result, while the crosses and circles show the NLO 
and N 2 LO results. The dotted lines indicate the pion masses at which a,^ 2 di- 
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verges because the second and third excited states of the triton appear at the 

1 /2 

neutron-deuteron threshold. These singularities in a ' . {M n ) are a clear signature 
that the limit cycle is approached in the critical region. 

A final answer on the question of whether an infrared limit cycle can be real- 
ized in QCD can only be given by solving QCD directly. In particular, it would 
be very interesting to know whether this can be achieved by appropriately tun- 



ing the quark masses in a Lattice QCD simulation 100 . The first full lattice 



QCD calculation of nucleon-nucleon scattering was reported in 97 but statisti- 
cal noise presented a serious challenge. A promising recent high-statistics study 
of three-baryon systems presented also initial results for a system with the quan- 
tum numbers of the triton such that lattice QCD calculations of three-nucleon 



systems are now within sight |101|. For a review of these activities, see Ref. 102 



Such calculations require a detailed understanding of the modification of the Efi- 
mov spectrum in a cubic box. For identical bosons, there are significant finite 



volume shifts even for moderate box sizes 103 . These shifts show universal scal- 



ing behavior which could be exploited to reduce the computational effort 104 
The extension of these studies to the triton case is in progress. 



3.3 Halo Nuclei 

A special class of nuclear systems exhibiting universal behavior are halo nuclei 



27, 105 . Halo nuclei consist of a tightly bound core surrounded by one or more 
loosely bound valence nucleons. The valence nucleons are characterized by a very 
low separation energy compared to those in the core. As a consequence, the radius 
of the halo nucleus is large compared to the radius of the core. A trivial example 
is the deuteron, which can be considered a two-body halo nucleus. The root 
mean square radius of the deuteron is about three times larger than the size of 
the constituent nucleons. Halo nuclei with two valence nucleons are particularly 
interesting examples of three-body systems. If none of the two-body subsystems 
are bound, they are called Borromean halo nuclei. This name is derived from 
the heraldic symbol of the Borromeo family of Italy, which consists of three rings 
interlocked in such way that if any one of the rings is removed the other two 
separate. The most carefully studied Borromean halo nuclei are 6 He and 11 Li, 
which have two weakly bound valence neutrons |l05| . In the case of 6 He, the 
core is a 4 He nucleus, which is also known as the a particle. The two-neutron 
separation energy for 6 He is about 1 MeV, small compared to the binding energy 
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of the a particle which is about 28 MeV. The neutron-a (na) system has no 
bound states and the 6 He nucleus is therefore Borromean. There is, however, 
a strong P-wave resonance in the J = 3/2 channel of na scattering which is 
sometimes referred to as 5 He. This resonance is responsible for the binding of 
6 He. Thus 6 He can be interpreted as a bound state of an o-particle and two 
neutrons, both of which are in P 3 / 2 configurations. 

Because of the separation of scales in halo nuclei, they can be described by 
extensions of the pionless EFT. One can assume the core to be structureless and 
treats the nucleus as a few-body system of the core and the valence nucleons. 
Corrections from the structure of the core appear in higher orders and can be 
included in perturbation theory. Cluster models of halo nuclei then appear as 
leading order approximations in this "halo EFT" . A new facet is the appearance 
of resonances as in the neutron-alpha system which leads to a more complicated 
singularity structure and renormalization compared to the few-nucleon system 



discussed above 1106 



The first application of effective field theory methods to halo nuclei was car- 
ried out in Refs. 106 107 , where the na system (" 5 He") was considered. It 



was found that for resonant P-wave interactions both the scattering length and 
effective range have to be resummed at leading order. At threshold, however, 
only one combination of coupling constants is fine-tuned and the EFT becomes 
perturbative. Because the na interaction is resonant in the P-wave and not in 
the S-wave, the binding mechanism of 6 He is not the Efimov effect. However, 
this nucleus can serve as a laboratory for studying the interplay of resonance 
structures in higher partial waves. 

Three-body halo nuclei composed of a core and two valence neutrons are of 
particular interest due to the possibility of these systems to display the Efimov 
effect |3j. Since the scattering length can not easily be varied in halo nuclei, 
one looks for Efimov scaling between different states of the same nucleus. Such 



analyses assume that the halo ground state is an Efimov state * 2 They have 
previously been carried out in cluster models and the renormalized zero-range 
model (See, e.g. Refs. [108 - IIP) ). A comprehensive study of S-wave halo nuclei 
in halo EFT was recently carried out in Ref. Ill . This work provided binding 
energy and structure calculations for various halo nuclei including error estimates. 



Confirming earlier results by Fedorov et al. 108 and Amorim et al. 109 , 20 C was 



#2 We note that it is also possible that only the excited state is an Efimov state while the 
ground state is more compact. This scenario can not be ruled out but is also less predictive. 
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Figure 11: The one- and two-body matter density form factors F c , F n , F nc , and 
F nn with leading order error bands for the ground state of 20 C as a function of 
the momentum transfer k 2 . 



found to be the only candidate nucleus for an excited Efimov state assuming the 
ground state is also an Efimov state. This nucleus consists of a 18 C core with spin 
and parity quantum numbers J p = + and two valence neutrons. The nucleus 
19 C is expected to have a ^ + state near threshold, implying a shallow neutron- 
core bound state and therefore a large neutron-core scattering length. The value 
of the 19 C energy, however, is not known well enough to make a definite statement 
about the appearance of an excited state in 20 C. An excited state with a binding 
energy of about 65 keV is marginally consistent with the current experimental 
information. 

The matter form factors and radii of halo nuclei can also be calculated in the 



halo EFT 111, 112 . As an example, we show the various one- and two-body 



matter density form factors F c , F n , F nc , and F nn with leading order error bands 
for the ground state of 20 C as a function the momentum transfer k 2 from 



111 



in Fig. 11 The theory breaks down for momentum transfers of the order of the 
pion-mass squared (k 2 ~ 0.5 fm -2 ). 

From the slope of the matter form factors one can extract the corresponding 
radii: 



F(k 2 ) 



1 



k 2 (r 2 ) + 



(44) 



Information on these radii has been extracted from experiment for some halo 



nuclei based on intensity interferometry and Dalitz plots [113] ■ Within the er- 
ror estimates, the extracted values are in good agreement with the theoretical 



predictions from halo EFT [ill] . For the possible 20 C excited state, the halo 
EFT at leading order predicts neutron-neutron and neutron-core radii of order 
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40 fm while the ground state radii are of order 2-3 fm. The theoretical errors are 
estimated to be of order 10%. Assuming a natural value for the effective range 
of the n- 18 C interaction, tq ~ 1/M W , next-to-leading order predictions for these 
radii have recently been obtained |114 . The leading order results were found to 
be stable under inclusion of the leading effective range corrections and the typical 
errors could be reduced to about 1-2%. 

Scattering observables offer a complementary window on Efimov physics in halo 
nuclei and some recent model studies have focused on this issue. In particular, in 



Refs. 115 116 the trajectory of the possible 20 C excited state was extended into 
the scattering region in order to find a resonance in n- 19 C scattering. A detailed 



117 



study of n- 19 C scattering near an Efimov state was carried out in 

The simplest strange halo nucleus is the hypertriton, a three-body bound state 
of a proton, neutron and the A. The total binding energy is only about 2.4 
MeV. The separation energy for the A, E\ = 0.13 MeV, is tiny compared to the 
binding energy = 2.22 MeV of the deuteron. The hypertriton can therefore 
also be considered a two-body halo nucleus. It has been studied in both two- 



body and three-body approaches 118-120 . A study of the hypertriton in the 



halo EFT was carried out in Ref. 121 . The AN scattering lengths are not 
well known experimentally since the few scattering data are at relatively high 
energies. If the AN scattering lengths are large, the hypertriton is likely bound 
due to the Efimov effect. In this case there will also be a correlation between 
the Ad scattering length a^d and the hypertriton binding energy analog to 



the Phillips line in the neutron-deuteron system 119 . In Fig. 12, we show this 
Phillips line correlation for three values of the AiV pole position 7^ |121|. For 
small hypertriton energies B^ , the different Phillips lines coincide exactly (the 
physical hypertriton corresponds to MB^ ~ 1.06 7^) and deviate from each other 
only at very large binding energies where the EFT breaks down. For all practical 
purposes the Phillips line is therefore independent of 7$. Whether the Efimov 
effect plays a role for the hypertriton is an open question. Most modern hyperon- 
nucleon potentials, however, favor a natural AN scattering length [52] . 

Another powerful method that can be used to investigate the Efimov effect 
in three-body halo nuclei at existing and future facilities with exotic beams 
(such as FAIR and FRIB) is Coulomb excitation. In these experiments a nu- 
clear beam scatters off the Coulomb field of a heavy nucleus. Such processes can 
populate excited states of the projectile which subsequently decay, leading to its 
"Coulomb dissociation" 1122] . The halo EFT offers a systematic framework for a 
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Figure 12: Phillips line in the hypertriton channel for different values of the AiV 
pole position 7; (all quantities are in units of the deuteron pole position 74 rj 45 
MeV). 

full quantum-mechanical treatment and can be used to predict the signature of 
the Efimov effect in these reactions. 



3.4 Three- Alpha System and Coulomb Interaction 

The excited + state in 12 C is known as the Hoyle state. Its properties are impor- 
tant for stellar astrophysics since it determines the ratio of carbon to oxygen in 
stellar helium burning. Efimov suggested that the Hoyle state could be explained 
as an Efimov state of a particles [3|[lT] . (For a more detailed discussion, see 
Ref. 1 123].) In this case, the universal properties arc modified by the long-range 
Coulomb interaction. The modified Efimov spectrum needs to be understood 
before any definite statement about the nature of the Hoyle state can be made. 

Many recent studies have focused on the consistent inclusion of the Coulomb 
interaction in two-body halo nuclei such as the pa and aa systems [124 125 



In particular, the aa system shows a surprising amount of fine-tuning between 
the strong and electromagnetic interaction. It can be understood in an expansion 
around the limit where, when electromagnetic interactions are turned off, the 8 Be 



ground state is exactly at threshold and exhibits conformal invariance 125]. In 
this scenario, the Hoyle state in 12 C would indeed appear as a remnant of an ex- 
cited Efimov state. In order to better understand the modification of the Efimov 
spectrum and limit cycles by long-range interactions such as the Coulomb inter- 
action, a one dimensional inverse square potential supplemented with a Coulomb 



interaction was investigated in 126 . The results indicate that the counterterm 
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required to renormalize the inverse square potential alone is sufficient to renor- 
malize the full problem. However, the breaking of the discrete scale invariance 
through the Coulomb interaction leads to a modified bound state spectrum. The 
shallow bound states are strongly influenced by the Coulomb interaction while the 
deep bound states are dominated by the inverse square potential. These results 
support the conjecture of the Hoyle state being an Efknov state of a particles 
but a full calculation of the 3a system including Coulomb in the halo EFT is 
missing. Calculations with the fermionic molecular dynamics model and electron 
scattering data, however, support a pronounced a cluster structure of the Hoyle 
state 127 



4 Applications in Particle Physics 



4.1 Hadronic Molecules 



In recent years many new and possibly exotic charmonium states have been ob- 
served at the B-factories at SLAC, at KEK in Japan, and at the CESR collider 



at Cornell. This has revived the field of charmonium spectroscopy 128-130 
Several of the new states exist very close to scattering thresholds, and can be 
interpreted as hadronic molecules. If they are sufficiently shallow, one may ask 
whether there are any three-body hadronic molecules bound by the Efimov effect. 
A particularly interesting example is the X(3872), discovered by the Belle 



collaboration 131 in — > K^ z tt + tt J/ip decays and quickly confirmed by CDF 
132 , DO 133 , and BaBar 134 . The state has likely quantum numbers J PC = 



1 ++ and is very close to the D*°D° threshold. As a consequence, the X(3872) 
has a resonant S-wave coupling to the D*°D° system. An extensive program 
provides predictions for its decay modes based on the assumption that it is a 
D*°D° molecule with even C-parity: 

(D*°D°) + = (D*°D° + D°D*°) . (45) 
v2 

This assumption naturally explains several puzzling features such the apparently 
different mass in the J/ipTr + ir~ and D*°D° decay channels and the isospin vio- 



lating decays [135 , 136 . A status report with references to the original literature 



can be found in Ref. 137 



Using the latest measurements in the J/^ir + ir channel, the mass of the 
A"(3872) is [138]: m x = (3871.55 ± 0.20) MeV, which corresponds to an energy 
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E x = (-0.26 ± 0.41) MeV . (46) 

The central value corresponds to a (D*°D°) + bound state with binding energy 
Bx = 0.26 MeV (but a virtual state cannot be excluded from the current data in 



the J/V-vr+vr- and D*°D° channels |139||i4l] ). The V(3872) is also very narrow, 
with a width smaller than 2.3 MeV. 

Because the X(3872) is so close to the D*°D° threshold, it has universal low- 



energy properties that depend only on its binding energy 142 . Close to thresh- 
old, the coupling to charged D mesons can be neglected because the D* + D~ 
threshold is about 8 MeV higher in energy. Therefore, the properties of the 
V(3872) can be described in a universal EFT with contact interactions only. Un- 



fortunately, there is no Efimov effect in this system 142 , so universal bound 
states of the X and D° or D*° mesons do not exist. The reason for this that 
there is not a sufficient number of pairs with resonant interactions as only the 
D*°D° and D°D*° interactions are resonant. However, it is possible to provide 
model-independent predictions for the scattering of D° and D*° mesons and their 
antiparticles off the X(3872). This scattering process is to leading order deter- 
mined by the D*°D° and D°D*° interactions only. 

The corresponding cross sections as a function of the center-of-mass momentum 



k obtained in |138j are shown in Fig. 13 The difference between the contribution 
of S-waves (L = 0) and the full cross section (including all partial waves up to 
L = 6) is negligible for momenta below the bound state pole momentum 7. Our 
results are given in units of the scattering length and may be scaled to physical 
units once a is determined. At present the error in the experimental value for 



Ex in Eq. (46) implies a large error in the scattering length. In particular, we 
obtain the ranges 7 = (0 . . . 36) MeV for the pole momentum and a = (5.5 . . . 00) 
fm for the scattering length with central values 7 = 22 MeV and a = 8.8 fm. 
Using the central value of the scattering length, we obtain for the scale factor 
a 2 = 0.78 barn. This factor can become infinite if the V(3872) is directly at 
threshold, while the lower bound from the error in Ex would give a value of 0.3 
barn. Even in this case the total cross section at threshold will be of the order 
300 barns for D°X scattering and 1000 barns for D°*X scattering. It may be 
possible to extract the scattering within the final state interactions of B c decays 
and/or other LHC events. Observation of enhanced final state interactions would 
provide an independent confirmation of the nature of the X(3872). 
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Figure 13: Total cross section for scattering of D° and D*° mesons off the 
X(3872) for S-waves (L = 0) and including higher partial waves with L < 7, in 
units of the scattering length a. The cross section is the same for the scattering 
of particles as it is for the scattering of antiparticles. 

There may also be hadronic three-body molecules that are bound due to the 
Efimov effect but currently no strong candidate states are known. This situation 
will be improved by new experiments at facilities such as FAIR and Belle II which 
have a dedicated program to study exotic charmonium states. 

5 Summary and Outlook 

Any few-body system with short-range interactions that has a two-body scat- 
tering length larger than the range of the underlying interaction will display 
universal properties and Efimov physics. This statement is independent of the 
typical length scale of the system and atomic, nuclear and particle physics can 
provide examples of universality. While much progress has been made recently in 
experiments with ultracold atoms, the concept of Efimov physics was originally 
devised for the few-nucleon problem. As we have shown, it can serve as a start- 
ing point for a description of very low-energy nuclear phenomena in an expansion 
around the unitary limit. 

In this review we have discussed the manifestation of Efimov physics with a 
strong emphasis on nuclear and particle physics. There are a number of such sys- 
tems that display low-energy universality associated with Efimov physics. The 
most important example from nuclear physics is the triton. The nucleon-nucleon 
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scattering length is large compared to the range of the internuclear interaction. 
Phase shift equivalent NN potentials will therefore necessarily give results for 
triton binding energy and neutron-deuteron scattering length that are correlated 
and lie on the Phillips line. The implications of universality on the four-nucleon 
system have been also been explored. It was found that the Tjon line, a corre- 
lation between three-nucleon and four-nucleon binding energies is a result of the 
large two-nucleon scattering length. 

Halo nuclei might provide a further example of few-body universality. In par- 
ticular two- neutron Halos such as 20 C could display Efimov physics in form of an 
excited three-body state due to the large core-neutron scattering length. While 
first studies of bound state observables have become available, scattering cal- 
culations represent an exciting opportunity for future applications of the EFT 
approach. In particular so-called p — t reactions in which a triton is formed in a 
collision of a proton and a two-neutron halo will provide an important benchmark. 

Several new charmonium states have recently been discovered close to scatter- 
ing thresholds and can be interpreted as hadronic molecules. If they are suffi- 
ciently shallow, they have universal properties associated with large scattering 
length physics. The best known example is the X(3872) which may be inter- 
preted as a D*°D° molecule with even C-parity. There may also be three-body 
hadronic molecules bound by the Efimov effect but currently no strong candidates 
are known. 

The separation of scales between scattering length and range facilitates the 
application of an EFT that reproduces at leading order the results obtained by 
Efimov and is known as the pionless EFT. Within this framework corrections 
to the zero-range limit can be calculated systematically in a small parameter 
expansion in powers of hi and £/a. The number of few-body calculations that 
include higher order corrections is growing. As a consequence, the expansion 
around the unitary limit provides a useful starting point for a controlled descrip- 
tion of very low energy phenomena in nuclear and particle physics. The intricate 
consequences of Efimov physics are explicit in this framework and universal cor- 
relations between observables arise naturally. Moreover, this theory is an ideal 
tool to unravel universal properties and establish connections between different 
fields of physics. 

The constituents of nuclear few-body systems can have charge in contrast to 
the neutral atoms used in experiments with ultracold gases. Electroweak observ- 
ables provide thus additional information on few-body universality. However, 
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they are also of interest by themselves and the pionless EFT guarantees a con- 
sistent framework for the calculation of observables with the minimal number 
of parameters as it is straightforward to consider external currents in any EFT. 
First calculations for observables such as form factors and capture rates have been 
performed but many more remain to be calculated. Of very high interest are in 
this context thermal capture rates in the four-nucleon sector that are relevant to 
big bang nucleosynthesis. 
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